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ACCELERATION  OF  AN  ELECTRON  RING  IN  A  MODIFIED 
BETATRON  WITH  TRANSVERSE  PRESSURE 

I.  Introduction 

Interest  in  the  modified  betatron  derives  in  part  from  its  promise  as  a 
high  current  accelerator.*-*  It  is  designed  to  overcome  the  current 
limitations  in  a  conventional  betatron,  where  the  defocusing  self-field  forces 
(which  scale  like  E/y  ,  where  E  is  the  self-electric  field,  and,  y  is  the 
relativistic  factor)  must  be  less  than  the  focusing  forces  from  the  external 
magnetic  field  index  n  (0  <  n  <  1).  The  conventional  betatron  then  has  two 
viable  Injection  regimes  -  one  where  the  injected  energy  is  low  and  the 
current  (or  density)  is  low,  and  the  other,  where  the  current  (or  density)  is 
high  and  the  injected  energy  is  high.  The  inclusion  of  an  external  toroidal 
field,  BQ,  in  the  modified  betatron  allows  injection  and  trapping  at  lower 
energy  and  high  current  (1-10  kA).  Other  theoretical  studies  of  the  behavior 
of  this  accelerator  have  either  used  particle  simulations  to  determine  the 
self-consistent  motion, or  have  assumed  a  given  density  and  current  profile 
and  calculated  the  motion  of  test  particles  in  fixed  fields.  *  *  *  The 
drawback  of  particle  simulations  is  the  large  computational  cost,  particularly 
because  the  time  scales  of  acceleration  and  self-flux  diffusion  are  on  the 
order  of  10“*  -  10“^  sec  while  the  time  scale  of  one  particle  cyclotron  orbit 
is  on  the  order  of  10-*0  sec.  The  analytic  technique  is  limited  because  of 
the  assumed  charge  and  current  profile. 

o.in 

In  our  previous  work,  a  theoretical  technique  is  developed  to  extend 

the  methods  mentioned  above  and  to  analyze  both  the  self-consistent 

equilibrium  in  the  modified  betatron  and  the  evolution  of  this  equilibrium  as 

external  parameters  are  varied  slowly  compared  to  a  drift  period.  These 
Manuacript  approved  Augutt  2,  1984. 


calculations  assume  a  cold  (zero  pressure)  beam.  In  this  paper,  we  expand  our 
theoretical  technique  and  analyze  the  evolution  of  equilibria  in  the  modified 
betatron  with  transverse  pressure. 

An  investigation  is  also  made  to  determine  if  the  presence  of  pressure 
causes  any  qualitative  changes  in  behavior  which  can  be  measured  easily  in  the 
experiment;  such  effects  may  indicate  a  useful  diagnostic  measurement  of  beam 
pressure.  The  effect  of  random  transverse  motion  appears  as  a  pressure 
gradient  term  in  the  force  balance  equation  in  our  formalism  in  a  limit  that 
is  applicable  to  current  designs  for  modified  betatrons. 

As  will  be  shown  in  the  section  on  our  new  formalism,  we  find  that  many 
of  our  earlier  results  still  hold  to  first  order,  namely  that  the  fluid 
canonical  angular  momentum,  Pg,  characterizes  the  equilibria  (because  many 
Important  fluid  quantities  are  functions  of  Pa,  and  because  surfaces  of 

constant  PQ  are  drift  surfaces  of  single  electrons).  The  pressure  H  manifests 

Itself  as  an  additional  defocusing  force  on  the  beam  (with  force  density 
-*  2II/a),  tending  to  increase  its  minor  radius  a.  The  pressure  also  produces 
an  outward  "hoop  force"  (with  force  density  ^n/R),  tending  to  Increase  the 
major  radius  R  of  the  beam.  These  two  effects  can  be  seen  in  our  results 
by:  a)  an  Increase  in  the  poloidal  x  Bt  drift  on  the  beam  (where  Fp 
represents  the  total  poloidal  defocusing  force,  and  Bt  the  toroidal  magnetic 
field),  and  b)  a  slight  Increase  In  the  vertical  field  B2  necessary  to  hold 

the  beam  at  its  equilibrium  major  radius  R. 

We  have  performed  two  different  classes  of  calculations.  The  first  class 
involves  a  series  of  separate  equilibria  with  similar  beam  energies  but 
various  levels  of  pressure.  The  other  class  involves  the  calculation  of  a 
series  of  adiabatically  coupled  equilibria  during  acceleration.  These  will  be 
discussed  in  Section  IV. 


In  Section  II  our  new  fluid  model,  including  pressure,  is  presented. 

This  is  followed  in  Section  III  with  a  description  of  how  the  new  model  is 
Incorporated  into  the  numerical  code,  as  well  as  a  qualitative  discussion  of 
the  effects  of  pressure  on  the  system.  Next,  we  discuss  our  results  in 
Section  IV  and  summarize  in  Section  V. 

II.  Fluid  Equations  and  Betatron  Ordering 

In  this  section  we  derive  fluid  equations  for  a  warm  relativistic 
electron  beam  in  a  modified  betatron.  We  derive  an  ordering  scheme  for 
diamagnetic  equilibria  which  indicates  the  conditions  under  which  the 
transverse  thermal  spread  of  the  beam  can  be  described  by  a  scalar  pressure, 
and  the  conditions  under  which  this  pressure  can  be  assumed  to  be  constant  on 
the  fluid  Pg  surfaces. 

The  first  moment  of  the  relativistic  Vlasov  equation  for  f(x,g,t)  is 

(n<5»  +  V«Jfvgd3p  -  nq  (E  +  i  <v>  x  B),  (1) 

which  can,  with  the  aid  of  the  continuity  equation  3n/3t  +  V»n  <v>  -  0,  be 
written  in  the  form 


n  <v>»Vg  ■  nq  (E  +  ■£•  <v>  x  B)  -  7»II, 

3  3 

for  equilibria  (3/3t  *0),  where  g  -  ymv,  n  ■  Jfd  p,  n<v>  -  J  fvd  p, 
n  <g>  -  jfgd3p, 


(2) 


m  <v> 

?T7^’  (3b) 


and 


nij  ’  id3p?sr  (pi  ■  V  (pi  '  V*  (4) 

A  manifestly  covariant  formulation  of  the  fluid  equations  has  recently  been 
described  by  Newcomb. ** 

The  class  of  equilibria  we  wish  to  describe  have  distribution  functions 
of  the  form 


f(H,P0,u)  -  W(PQ,y)  6[H  -  E(P0) J ,  (5) 

2 

where  H  is  the  particle  energy  ymc  +  q<j»,  PQ  -  ymrvQ  +  qi|»/c  is  the  canonical 
angular  momentum  (41  ■  rA0),  and 

<Ei  ■  ip2 

“‘S-e’  '  -2B9(rTzV  <6> 

is  the  first  adiabatic  invariant.  Here,  (r,0,z)  is  a  cylindrical  coordinate 
system  with  z  the  vertical  direction  and  9  referring  to  the  toroidal  symmetry 
direction.  Also,  1  refers  to  components  perpendicular  to  the  toroidal 

A 

magnetic  field  B00.  The  conditions  under  which  y  is  invariant  are  discussed 
in  Appendix  A.  This  distribution  function  has  a  distribution  of  magnetic 
moment  y  and  contains  the  cold  beam  case  of  Refs.  8  and  9  as  a  special  case 
W(P0,y)  -  Wo(P0)6(y). 

In  general,  a  distribution  such  as  (5)  cannot  be  required  to  have  only 
transverse  thermal  spread,  in  spite  of  the  delta  function  in  energy,  because 


4 


the  distribution  of  y.  That  is,  no  element  of  the  stress  tensor  g  Riven  in 
Eq.  (4)  is  exactly  zero  unless  the  distribution  is  cold.  Nevertheless,  there 
is  an  ordering  scheme  appropriate  to  the  modified  betatron  in  which  g  is 


isotropic  in  the  poloidal  plane  (r,z).  Defining  u)p  to  be  the  plasma  frequency 

and  at  to  be  the  (nonrelatlvlstlc)  cyclotron  frequency  in  the  toroidal  field, 
c 

the  externally  imposed  vertical  field,  and  ut  *  eBv/mc,  we  take 


2 

!!e 

2 


u 

c 


0(1), 


a 

IT 


0(e2), 


(7a) 

(7b) 


Y  -  0(e"1-6), 


0(e1"6), 


0(1), 


jP 

a 


0(e2/3). 


(7c) 

(7d) 

(7e) 

(7f) 


Here,  p  is  the  electron  gyroradius  in  the  toroidal  field,  a  is  the  minor 

radius  of  the  beam  and  R  the  major  radius.  Since  this  ordering  must  describe 

the  beam  from  Injection  (y  -  4)  up  to  the  diamagnetic-paramagnetic 

2 

transition  (y  >  15),  and  since  the  ratio  of  eE/y  to  focusing  forces  scales 
3 

as  1/y  ,  the  extra  parameter  6  is  required  for  the  ordering  to  be  valid 
throughout  the  acceleration.  For  parameters  relevant  to  the  NRL  modified 
betatron,  e  **■  1/10;  also,  S  -  -1/3  corresponds  to  injection  and  5  -  1/3 
corresponds  to  transition. 


Note  that  the  scalings  with  5  in  Eq.  (7)  reflect  the  invariants  of  motion 

during  acceleration  or  any  adiabatic  evolution  of  the  beam.  In  particular, 

2  2 

<i)  /(i>  is  Independent  of  time  (i.e.  of  5)  because  conservation  of  the  total 
P  c 

Q  Q 

toroidal  flux  #  through  the  beam  »  implies  that  the  density  remains 

approximately  constant,  and  the  acceleration  scheme  we  are  concerned  with  does 

not  Involve  changing  the  external  toroidal  magnetic  field.  (The  self  magnetic 

field  is  smaller  in  this  ordering.)  Similarly,  relation  Eq.  (7e)  represents 

the  fact  that  y  is  proportional  to  By  during  acceleration,  due  to  PQ 

conservation.  Finally,  since  p  is  the  flux  through  the  gyro-orbit,  relation 

Eq.  (7f)  follows  from  p/a  **  (p/$)  .  Also,  Eq.  (7f)  guarantees  that  the 

electron  velocity  about  the  guiding  center  put  fy  is  of  the  same  order  as  the 

c 

2 

drift  velocity  cE/y  B  at  injection,  a  plausible  assumption  for  most  injection 
schemes.  In  fact,  for  this  ordering  we  find  that  the  drift  and  gyration 
velocities  obey  v^/c  -  and  v^/c  *  With  this  ordering,  we  find 

that  the  forces  eE/y  ,  ev^Bg/c,  the  external  focusing  forces,  and  the 
pressure  gradient  are  in  relation 

el+3«:  el+36:  £2:  e2*4/3.  (*} 

We  see  that  the  first  two  forces  dominate  at  injection  (5  *  -  1/3),  whereas 

all  four  are  of  the  same  order  at  transition  (6  -  1/3).  Note  that  for 

intermediate  times,  -1/3  <5<l/3,  the  pressure  gradient  dominates  focusing 

2  1+6 

forces.  Also  note  that  q$/ymc  v/y  *«•  1/y  -  e 

Using  this  ordering  we  first  estimate  the  distribution  of  y  which  is 
due  to  the  distribution  of  u  shown  in  Eq.  (5).  Assuming  the  energy 
E(P0)  -  Eq  +  n(PQ  -  P0q)  as  in  Ref.  8,  H  -  E  (P0)  from  Eq.  (5)  and  the 
definition  of  y  give 


Yoc  -  YmrOv0  -  Hfl(r,z)  =  Eq  -  0PQo  +  qfl'F/c  -  q$. 


2  i  a.  2/  2  2  ^  ^  n2,  2  2 

Y  -  1  +  p9/m  c  +  (gid  +  Si  )  /m  c  , 


where  p0  ■  Ymv0,  5j_<j  1®  c^e  momentum  associated  with  the  guiding  center  drift, 

and  g  is  the  momentum  associated  with  gyration  about  the  guiding  center. 

Prom  Eq.  (7)  we  find  p^  **  e  ,  p*  e  '  ,  the  latter  constant  since  it  is 

essentially  magnetic  moment.  Putting  Eqs.  (9)  in  dimensionless  form  by 

2 

defining  q0  -  p0/mc,  ^  -  g^/mc,  K  -  HQ(r,z)/mc  ,  and  u>  -  rO/c  and  combining, 
we  find 


1  +  qe  +  qi  "  <K  +  “V2, 


Expanding  for  q^«l,  K»1  [from  Eq.  (7),  these  quantities  scale  as 

6*^  and  e  *  respectively],  i.e.  treating  1  +  qj  as  a  perturbation  in  Eq. 

(10),  we  find 


9o  1  -  u>’ 


(11a) 


1  +  qI 

6q9  "  2K  ’ 


(lib) 


die) 


«Y  -  - 


<*>(1  +  qp 


(lid) 


Therefore,  the  thermal  contribution  to  5q  and  5y,  5qa.  *  ~q,  / 2K  and 

9  “  t  Ig 


(12) 


6y.  -  -  a*i ,  / 2K,  satisfy 
t  ig 


«0o  \ 


-  e 


10/3  +  25 


Also 


5v0  6<l0 


q9o5Y 


(1  -o»)2(l  +  qj) 
2K5 


(13) 


the  thermal  contribution  to  which  has 


5v 


0t 


c 


10/3  +  26 
&  • 


(14) 


2  2  25  2  a  /o 

Using  -  e  ,  q*g  e  7  ,  we  find  that  the  drift  contributions  to 

6qQ,  6y  and  6v0  are  smaller  than  the  thermal  contributions  for  S>  -1/3.  The 
issue  of  conservation  of  y  in  this  ordering  is  addressed  in  Appendix  A. 

We  now  use  the  relations  Eqs.  (11)-(14)  to  evaluate  the  stress  tensor  n 
of  Eq.  (4)  for  equilibrium  distribution  functions  of  the  form  Eq.  (5).  First, 
noting  d  p  »  (Bg/r )  dy  dx  dPQ,  where  x  “  tan  [(pz  -  P2>/(Pr  -  Pr>J  la  the 
gyrophase,  we  find  that  the  density  is  given  by 


n 


r 


jdydx  } 


dP9W(P0,y) 


5  [P0  -  P0(r,r,p2)J 
v0/r  --i?(P0) 


(15) 


A 

where  P0  »  ?g(r,*,pp  is  the  solution  of  H  ■  E(P0),  found  in  Eqs.  (9)  -  (11) 
for  the  special  case  ft  ■  dE/dPQ  ■  const.  For  this  case,  we  find 


(16) 


P0(r,z,p2)  -  P0°^(r, z)  +  mcr5q0, 

where  P^  neglects  the  drift  and  thermal  corrections  in  Eq.  (9).  From  the 
above  ordering,  the  thermal  and  drift  contributions  to  PQ  and  v0  are 
negligible  and  we  find 

2*B  •  /«\ 

n  - J  '(r,*),  uj-  (17) 

0o  o 

Also,  using  this  ordering  we  find  that  g  *  <p>  *  yq»<v>  in  Eq.  (12)  and  that 
the  two  diagonal  poloidal  components  of  the  stress  tensor  are 


n  -  ji 

rr  zz 


2  it  b; 


*  r~.fr-  rST  I  ydyW[Pg0>(r, *),nj , 

O  So  o 


(18) 


while  the  off  diagonal  terms  are 


2Be  dpe 

-  ~  Jdudx  8ln  x  cos  x  ^y(v  -ray  w(Vw) 

9 

X  5[P0  -  P0(r,z,p2)j.  (19) 


If  we  neglect  the  thermal  and  drift  corrections  to  P0,  y»  and  v0,  these  terms 
vanish  by  the  integration  over  the  gyrophase  x*  similar  considerations  we 
find  nQr  -  n0z  -  nr0  -  nz0  -  n00  -  o.  Thus,  -  vn(r,z),  with  n  -  nrr 


9 


Using  these  conclusions,  we  find  that  the  toroidal  component  of  the 


equation  of  motion,  Eq.  (2) 


m0  . 

nv*  Vp„  +  -  p 

■*  K9  r 


— -  n  — -  n  -  —  i 

rc  '  w  3r  ur9  3z  z9  r  1 


r9* 


(20) 


gives 


Ip*  7^omr70o  +  “  °* 


(21) 


Since  the  poloidal  velocity  v^  is  }^/ne  “  (c/4irne)7g  x  70 
(B  ■  7\|)  x  79  +  g70) ,  and  since  the  quantity  in  parentheses  in  Eq.  (21)  is 
we  find,  as  in  the  cold  fluid  case.  Ref.  9 

0 


8 


8 


(22) 


The  poloidal  components  of  the  equation  of  motion,  neglecting  poloidal 
centrifugal  force,  are 


2  *,B(0). 

Y  mv.  qv  gg  (?'  ) 

0  -  -2 — §2.  7r  -  q7<J>  +  — -  ^ - ® —  7T>£0) 

r  rc  W2  9 


7n 

n 


(23) 


-  -r  <’  -  "*9o,<w.>  -  -  -S£-t7p9°>  - 

4xnr 


To  lowest  order  in  e,  II,  given  by  Eq.  (18)  is  a  function  of  Pa^. 

2  2  2  2 
Specifically,  the  relative  variation  of  B  -  g(P  )  /r  scales  as  5r/r  -  e  , 

0  0  o 

2 

and  the  relative  variation  of  y  is  eA$/ymc  ,  where  A$  is  the  variation 

of  $  along  a  PQ  surface.  For  highly  elongated  drift  surfaces,  as  can  occur 
Q  Q  —1  14-6 

j  near  transition,’  the  quantity  scales  as  y  -  e  {See  Appendix  A).  For 

I 

circular  flux  surfaces,  the  relative  variation  of  $  along  a  P.  surface  (due  to 

0 


the  toroidally  generated  displacement  between  these  surfaces)  is  of  order 
2  3+5 

a/R  "  e  ,  so  6y/y  scales  as  e  .  Finally,  the  denominator  v,  -rQ  scales 

7  O 

2 

as  Sr/rQ  **  e  as  long  as  the  two  terms  which  make  up  the  denominator  do  not 

nearly  cancel.  Thus  the  relative  change  in  pressure  along  a  P  surface  is  at 

9 

1+6 

worst  e  .  In  fact,  away  from  transition,  the  flux  surfaces  are  nearly 

always  circular  so  the  relative  pressure  variation  is  actually  much  less. 

Since  the  pressure  Itself  is  not  a  dominant  effect  except  near  transition 

(see  Eq.  (3)),  the  variation  of  II  along  a  P0  surface  is  negligible.  Thus 

2 

n(r,z)  -  n(PQ).  Using  these  relations,  which  imply  v.  7(vA  y  )  c  Vy  and 
0  0O  00  o  -  o 

VII  ■  H"(Pg0^)7Pg^\  we  find  (dropping  the  o  superscripts) 


-  £L_j  vp 
dP„J  e 


M'(P0) 

4nnr 


n'(pQ) 


J 


(24) 


We  can  use  this  equation  to  solve  for  the  density,  obtaining 


n 


vQ  -  rJ2(PQ) 


gg'(Pfl) 

l— — t~  +  n'(P0)J. 


4wr 


(25) 


Note  that  this  reduces  to  the  non-relativistic  MHD  result 
2 

j0  ■  cr [gg'(ij>)/4irr  +  H'(i|>)J,  obtained  from  ^  x  B  »  cVJI,  when  inertia  and 
electric  field  vanish  (P-+  q\p/c,  fl  ♦  0). 

It  is  worth  noting  that  H  can  be  specified  to  be  an  arbitrary  function  of 
P0.  Indeed,  using  the  approximations  under  which  II  is  a  function  of  P0  alone 
we  can  write 


n(PQ)  -  Y(P  )  J  ydpW(P9,  u  ), 
o 


(26) 


2 

where  Y(P0)  is  a  drift  surface  average  of  2uB9/[yom(v9o-  rQ)j.  We  can  find  a 


large  class  of  functions  W(Pa,  u)  that  produce  n(Pa),  for  example 

0  0 

W(P9,w)  -  WoCP0>5[vi  -  U(?Q)J,  (27) 

with  the  condition  “  n(P0)/Y(P0). 

I 

III.  Incorporation  of  Formalism  Into  Numerical  Code 

In  this  section  we  will  discuss  the  Incorporation  of  our  new  pressure 

formalism  Into  the  original  computer  code.  The  calculation  of  a  single 

equilibrium  (or  the  first  In  a  series  of  equilibria  adlabatlcally  connected) 

R— in 

proceeds  as  In  our  previous  work,  A  except  that  the  density  Is  calculated 
using  Eq.  (25).  For  convenience  we  model  the  Initial  density,  pressure 
gradient,  and  K(PQ)  ■  gg>(Pg)  as  linear  functions  of  PQ; 

Wty  ■  W  '  V/<p8c  -  %.>*  (2S) 


and 


n'(pa)  -  -2n  (pn  -  pa)/(ptt  -  pn  )  , 

v  0  ov  0c  0  0c  0o  ’ 


(29) 


where  P0c  Is  the  value  of  P0  at  the  edge  of  the  beam,  and  P0q  Is  the 

(conserved)  minimum  value  of  Pa  at  the  center  of  the  beam.  Equation  (29) 

0 


gives  a  quadratic  pressure  profile  In  Pg, 


n(P0)  -  no 


<P0c  -  V 
<P0c  -  P0o> 


2  * 


(30) 


Since  the  maximum  density  n(Pa  )  and  maxlraun  pressure  II  are  determined  from 


the  input  parameters  (of  injected  current  and  peak  pressure),  the  constant  KQ 
can  be  determined  from  Eqs.  (25),  (28)-(29).  (Note  that  as  in  our  previous 
work,  we  choose  fl(P0)  ”  constant.) 

By  integrating  Eq.  (25)  over  drift  surfaces,  we  obtain 


dN 

dpe 


1  ,  dl  1 

I88  qwjr 


2irIT(P0) 


(31) 


where  the  quantity  (dN/dP0)SP0  is  the  (conserved)  number  of  particles  between 
drift  surfaces  labeled  by  P0  and  P0  +  6P0.  In  our  iterative  solution  of  the 
equilibrium  equations  defining  the  system  we  conserve  the  initial  values  of 
dN/dP0,  the  beam's  equilibrium  major  radius  R  (location  of  VP0  -  0),  the 
minimum  P0Q,  and  the  toroidal  flux  in  the  beam  #t.  During  the  adiabatic 
evolution  of  these  equilibria,  the  change  in  pressure  must  be  determined  as 
external  parameters  change  adlabatically. 

We  begin  with  Eq.  (18),  defining  the  pressure  in  our  ordering  scheme 


To  lowest  order,  H  -  II(P_ ) ,  but  since  B„,  and  (v„  -  rfl)  vary  along  a  P- 

O  0  j  O  0 

surface,  let  us  redefine  the  pressure  as 

2itB9 

n<V  -  Vm(v-:  -rny  >  *i  W<V^dw’  (32> 

o  00  r0  o 

where  <  >„  is  an  average  over  a  PQ  surface. 

P9  9 

The  integral  in  Eq.  (32)  does  not  change  in  time  (p,  and  P  are  constants 
of  the  motion),  while  the  term  inside  the  brackets  changes  as  yq*  8,  SI,  and 


the  shape  of  P.  surfaces  change  in  time.  Using  Eq.  (32)  we  can  determine  the 


evolution  of  1  (?0,t)  in  time  from 


2  <  - T1 - > 

n(P0,t)  8  (P0,t)  Vt)r  (00o(t)  -  Ul(t))  P0(t) 

li?v3T‘7^  <  : nz . =7 


(33) 


YQ(0)r  ($0<>(O)  -  .,,(0))  p0(O) 


where  g  ■  rB0,  and  to  ■  rfl/c.  During  acceleration,  gg'(Pg)  is  determined  from 
Eq.  (31),  g  is  found  by  integration  of  gg",  and  Eq.  (33)  gives  1  (P0  ,t) . 

Since  gg'  itself  depends  on  1  (P0)  in  Eq.  (31),  the  cycle  must  be  Iterated 
until  it  converges. 

An  interesting  consequence  of  Eq.  (33)  is  that  pressure  actually 

Increases  during  acceleration.  A  simple  theory  for  an  infinitely  long 

cylindrical  beam  would  specify  p^  as  a  constant  of  motion  during  acceleration 

by  a  longitudinal  electric  field.  Therefore  it  would  predict  n  -  p^v^ 

2 

“P^/T  "v/y,  or  H  '  l/y.  It  can  be  seen  from  Eq.  (33)  however,  that  p  actually 
scales  as 


n '  Set)  (6e(t)  -  «(t))  y?Qm 


(34) 


During  acceleration  y  increases,  but  (0Q(t)  -  u>(t))  decreases  faster,  with  the 

C7 

net  effect  that  the  pressure  increases.  This  Increase  in  pressure  arises  from 
the  fact  that  the  parallel  and  perpendicular  parts  of  the  momentum  are  coupled 
by  the  pololdal  field  and  toroidal  geometry.  The  tendency  of  00  -  u>  to 
decrease,  i.e.  for  £2  =  dE/dPQ  to  increase,  is  due  to  toroidal  flux  invariance 
as  discussed  in  Ref.  9.  It  will  be  seen  in  our  acceleration  results  that  the 
Increasing  pressure  drives  the  pololdal  drift  of  the  beam  toward  the  Brlllouin 
limit  as  the  beam  accelerates.  A  brief  explanation  of  how  the  pressure 
affects  this  drift  may  be  useful. 


/.  \ 
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Let  p  denote  distance  from  beam  center.  Then  assuming  a  quadratic 


dependence  of  Pg  on  p  near  the  reference  orbit  we  can  write 


pe  "  p9o  +  (pec  “  peo)l>2/‘2- 


(35) 


where  a  is  the  beam  radius. 

Substituting  Eq.  (35)  into  Eq.  (30)  gives 


n  (p)  -  H0(1  -  p2/a2)2  . 


(36) 


The  pressure  gradient  plays  a  defocusing  role  in  poloidal  force  balance  as  can 
be  seen  by  writing  the  poloidal  fluid  equation  of  motion  near  the  reference 
orbit 


na 


(3  7) 


2 

where  qE/y  represents  the  sum  of  the  self  electric  and  self  magnetic 

q  q 

(Vg  x  Bge^)  forces  *  and  v  is  the  poloidal  velocity.  Equation  (3  7)  can  be 
rewritten  In  terms  of  the  poloidal  rotational  frequency  «  v/p  as 

q  2  2  ®  T  K  n 

2Y  %  +  +  — y-2  -  2T  n>cnD  -  °,  (38) 

p  na  m 


where  is  the  plasma  frequency,  and  wc  is  the  (non-relativistic)  cyclotron 
frequency  in  the  toroidal  field.  From  the  solution  of  Eq.  (38) 


it  can  be  seen  Chat  Che  maximum  (real)  value  of  ft-  is  still  ft_  ■  uj  /2-y ,  the 

u  DC 

Larmor  frequency  at  Brlllouin  flow,  buc  Chat  Che  pressure  serves  to  Increase 

2 

ftp  in  the  slow  mode,  effectively  adding  diamagnetic  drift  to  the  E  x  B/y 
drift. 

IV.  Results 

A.  Decoupled  Pressure  Equilibria 

In  this  first  set  of  calculations,  we  hold  certain  macroscopic 
quantities  fixed  and  find  beam  equilibria  for  different  values  of  pressure. 

At  each  level  of  pressure,  the  beam  energy  at  its  center,  the  total  number  of 
particles,  the  peak  density,  external  toroidal  field,  and  the  toroidal  flux  in 
the  beam  are  the  same,  but  the  calculations  are  not  adiabatically  coupled. 

Two  sets  of  calculations  are  discussed  here.  In  one  case  a  1  kA  beam 
with  y  ■  4,  minor  radius  a  -  2  cm,  and  major  radius  R  ■  100  cm  is  Injected 
into  a  toroidal  conducting  chamber  with  circular  cross  section  of  minor  radius 
12.5  cm,  and  major  radius  100  cm  (see  Fig.  1).  The  applied  vertical  field  Bz 
is  initially  set  equal  to  71.46  gauss  but  adjusts  Itself  slightly  to  maintain 
the  beam  at  its  major  radius.  We  set  the  field  index  n  ■  0.5  at  beam 

center.  The  applied  toroidal  field  B  is  held  at  2000  gauss.  We  varied  the 

0 

pressure  from  80  to  320  dynes/cm  .  As  a  point  of  reference  we  estimate  a 

pressure  which  corresponds  to  the  Lawson-Penner  emittance  e^_p  “  0 .3/1  (kA) 

2 

rad-cm.  Writing  pressure  in  terms  of  the  transverse  velocity,  H  *  fmnv^,  and 

relating  it  to  the  normalized  "emittance"  ■  ySa  -v ^ /v q ,  (vQ  ■  beam  velocity, 

2  2  2  2 

a  -  beam  radius)  we  get  H  “  nmv  e  ty 8  a  .  Setting  e  ■  eT  _  gives  nT  n  * 

on  n  L*r  u  ' 

187  dynes/cm  for  the  1  kA  beam.  In  an  urmagnetized  system,  this  emittance 
corresponds  roughly  to  random  transverse  energy.  In  a  magnetized  system  with 
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an  immersed  cathode,  this  random  transverse  energy  Is  manifested  as  finite 

Larmor  radius  (non-zero  p)  for  the  electrons. 

The  only  quantity  which  responded  significantly  to  this  pressure 

variation  is  the  poloidal  drift  averaged  across  the  beam.  This  term  is 

plotted  in  Fig.  2(a)  as  a  function  of  pressure.  In  Fig.  2(b)  we  show  the 

average  Larmor  radius  associated  with  the  pressure  (this  quantity  is  derived 

2 

from  the  simple  relations  n  ■  ymnv .  and  r.  ■  -ymcv./eB.).  Note  that  the 

1  D  10 

poloidal  drift  of  the  beam  is  calculated  from  V  x  B  «  4xJ/c,  giving  a  drift 

velocity  Vj  ■  c(Vg  x  9)/(4iraqr)  in  our  notation.  The  drift  frequency 

associated  with  v,  is  averaged  across  the  beam  and  normalized  to  the  Larmor 
a 

9  —1 

rotational  frequency  at  the  Brillouin  limit  ■  4.41  x  10  sec  . 

Other  quantities  such  as  the  vertical  magnetic  field  Bz  necessary  to 
maintain  R,  the  toroidal  field  ^  at  beam  center  the  poloidal  flux  at  beam 
center,  and  the  profiles  of  the  electrostatic  potential  $  and  poloidal 
flux  i|>,  changed  less  than  0.2%  over  this  range  of  emlttances.  The  vertical 
field  Increased  to  compensate  for  the  increased  hoop  forces,  and  the  BQ  field 
decreased  because  of  the  faster  diamagnetic  drift. 

In  our  second  case,  we  use  a  10  kA  beam  with  the  same  geometrical 
features  (a,  R,  etc.).  The  vertical  field  is  initially  120.1  gauss  and  the 
applied  toroidal  field  is  set  again  at  2000  gauss.  This  time  we  varied  the 
pressure  from  1.  x  10^  to  9.  x  10^  dynes/cm2,  where  -  21.5  x  10^ 

dynes/cm .  The  Larmor  frequency  in  this  10  kA  case  is  4.5  x  10  7  sec  ;  it  is 
slightly  larger  than  the  1  kA  case  because  of  a  space  charge  depression  of  y. 

As  in  the  previous  case,  the  poloidal  drift  Increases  with  the 
pressure.  The  drift  frequency  is  plotted  in  Fig.  3(a),  with  the  Larmor  radius 
in  Fig  3(b).  All  of  the  other  diagnostic  quantities  mentioned  above  vary  less 


The  quantity  which  depends  most  sensitively  on  transverse  pressure  is 
the  poloidal  rotation  frequency  of  the  electron  cloud.  This  indicates  that 
measuring  this  rotation  frequency  could  be  a  very  useful  diagnostic  of 
transverse  pressure.  In  HIPAC  the  rotation  frequency  was  measured  and  used 
as  a  diagnostic  for  total  charge.  In  the  modified  betatron,  because  of  the 
cancellation  of  self  electric  y  magnetic  field,  the  rotation  frequency  is  much 
less  sensitive  to  charge  and  much  more  sensitive  to  pressure.  Thus  if  the 
total  charge  could  be  measured  in  a  different  way  (say  by  laser  scattering), 
then  the  rotation  frequency  could  be  used  as  a  diagnostic  for  pressure. 

B.  Adiabatic  Acceleration  with  Pressure 

In  this  next  series  of  calculations,  we  examine  the  effect  of 
pressure  on  beam  acceleration  in  a  modified  betatron.  The  evolution  of  the 
pressure  during  acceleration  is  calculated  self-consistently  according  to  the 
formalism  described  In  Section  II.  In  other  respects,  however,  the  techniques 
used  here  follow  those  of  Refs.  8-10. 

We  compare  several  sets  of  calculations  emcompassing  a  wide  range  of 
injection  currents  (1-10  kA)  and  pressures  (0  -  "Lawson-Penner").  As  before, 
we  inject  a  y  ■  4  beam  with  minor  radius  a  ■  2  cm  and  major  radius  R  ■  100  cm 
into  a  toroidal  conducting  chamber  with  minor  radius  12.5  cm  and  major  radius 
100  cm.  A  variable  external  vertical  Bz  field  is  used  to  accelerate  the  beam 
and  hold  it  at  its  equilibrium  radius,  while  an  external  toroidal  B^  field  (of 
2000  gauss)  is  used  to  provide  beam  equilibrium  during  the  early,  diamagnetic 
stage  of  acceleration.  Our  calculations  follow  the  acceleration  until  the 
edge  of  beam  comes  close  to  reversing  its  diamagnetic  rotation. 

In  each  case,  we  choose  the  initial  pressure  so  that  the  beam's 
rotational  frequency  remains  below  the  Brlllouin  limit  (the  Larmor  frequency) 


throughout  the  acceleration  cycle.  This  Is  to  ensure  that  our  underlying 

assumption  of  small  pololdal  centrifugal  forces  remains  valid. 

In  Fig.  4,  we  plot  the  rotational  frequency  of  the  beam 

(calculated  by  averaging  across  the  beam)  during  acceleration  for  various 

injection  currents.  The  beams  were  injected  with  I  *  1  kA,  3  kA  and  10  kA, 

and  with  Initial  pressures  of  100,  800  and  8000  dynes/cm2  respectively.  For 

•1 

reference,  the  "Lawson-Penner  pressure"  for  each  case  is  187,  1.80  x  10  ,  and 
21.5  x  10-*  dynes/cm2  respectively.  (These  do  not  scale  exactly  as  I2  because 
the  peak  density  obtained  from  the  code  does  not  scale  exactly  as  I.)  In 
addition,  the  peak  pressure  as  a  function  of  beam  energy  Is  plotted  in  Fig.  5 
for  these  Injection  currents.  The  rotational  frequencies  of  Fig.  4  are 
plotted  as  ratios  of  the  Brlllouin  limit  frequency  8-  ■  eB./(2mYc)  which 
decreases  as  the  beam  accelerates,  but  the  trend  towards  the  Brlllouin  limit 
is  a  characteristic  of  beams  with  pressure.  Zero  pressure  beams  rotate  with 
frequencies  that  tend  away  from  the  Larmor  frequency  as  the  beam  accelerates 
(see  Fig.  6).  The  Increase  in  peak  pressure  with  acceleration  Is  responsible 
for  these  effects.  Higher  pressure  gradients  (near  the  center  of  the  beam) 
produce  larger  defocuslng  forces  and  hence  a  faster  drift  (?  x  S  ) 

p  9 

in  the  Bg  field.  As  discussed  in  Section  III,  the  Increase  of  pressure  with 
acceleration  is  an  unexpected  result  with  a  rather  subtle  cause  (see 
Eq.  (34)). 

The  Larmor  radius  of  the  beam  is  not  calculated  independently  In  our 

2 

formalism,  but  from  the  conservation  of  magnetic  moment  u  ■  irr  B  we  conclude 

U  0 

that  r^  changes  very  little  during  acceleration  because  the  relative  change 
in  Bg  dut  to  pololdal  beam  current  Is  small. 

One  interesting  result  we  have  observed  Is  that  the  diamagnetic 
transition  appears  to  occur  at  about  the  same  energy  for  beams  with  pressure 


as  those  without,  despite  the  larger  average  pololdal  drift  for  “warm" 

beams.  For  a  10  kA  beam  injected  at  y  ■  4,  for  instance,  we  find  by 

extrapolation  that  the  edge  of  the  beam  undergoes  the  transition  at  about 

y  ■  13.0  for  both  the  warm  [H  (t  ■  0)  ■  8  x  10  dynes/cm  I  and  the  cold  beam. 

1  o 

There  are  several  reasons  for  this  apparent  anomaly.  First  of  all,  we  note 

that  the  transition  from  diamagnetic  to  paramagnetic  rotation  is  determined  in 

our  calculations  by  following  the  separatrlx  of  the  PQ  surfaces  as  it 

approaches  the  beam  from  the  region  outside  the  beam.  In  this  region,  the  P0 

surfaces  are  also  outside  the  direct  influence  of  the  pressure  since  pressure 

vanishes  outside  the  beam.  Indirectly,  the  presence  of  pressure  Influences 

the  region  outside  the  beam  because  pressure  produces  an  outward  hoop  force 

which  increases  the  vertical  field  necessary  to  hold  the  beam  at  R,  which  in 

turn  increases  the  index  focusing  forces  throughout  the  system.  This  effect 

is  very  small  as  noted  in  our  decoupled  pressure  section.  Secondly,  note  that 

the  effect  of  pressure  within  the  beam  depends  on  distance  p  from  beam 

center.  As  seen  in  Eq.  (36),  our  pressure  has  a  quart ic  profile 
2  2  2 

n«IIo(l-p/a)  .  With  this  profile  the  pressure  gradient  is  largest 

inside  the  beam  and  drops  off  to  zero  near  the  edge,  producing  larger 

diamagnetic  drift  frequencies  at  the  center  of  the  beam  than  near  the  edge. 

According  to  these  arguments,  and  because  the  defocusing  effect  of 

2 

self-electric  and  self-magnetic  fields  (■»  E/y  )  is  essentially  the  same  for 
beams  with  and  without  pressure,  we  see  that  the  diamagnetic-paramagnetic 
transition  can  occur  at  the  beam's  edge  at  roughly  the  same  beam  energy. 
However,  the  pressure  will  undoubtably  affect  the  actual  drift  reversal  of  the 
particles  inside  the  beam,  which  occurs  after  transition. 


V. 


Conclusions 


In  our  analysis  of  pressure  effects  on  the  modified  betatron  we  have 
assumed  that:  a)  the  poloidal  drift  in  the  toroidal  field  is  below  the 
Brillouin  limit,  b)  that  the  parameters  affecting  the  beam  change 
adiabatically,  and  c)  that  the  parameters  of  the  experiment  obey  the  ordering 
in  Eq.  (7).  We  find  that,  as  expected,  the  pressure  increases  the  poloidal 
drift  frequency  of  the  beam  because  of  the  larger  defocusing  force.  In  spite 
of  this  increase,  however,  the  diamagnetic  to  paramagnetic  transition  begins 
at  about  the  same  energy  as  a  cold  beam  during  acceleration  because  of  the 
limited  scope  of  the  pressure  outside  the  beam.  Although  our  calculations 
cannot  treat  the  diamagnetic  to  paramagnetic  transition,  we  speculate  that  the 
actual  transition  process  will  be  much  different  in  the  warm  beam  case.  We 
also  find,  contrary  to  simple  estimates,  that  the  pressure  Increases  rather 
than  decreases  during  acceleration. 

Finally,  we  remark  that  the  sensitive  dependence  of  the  poloidal  drift 
frequency  on  the  pressure  could  be  used  as  a  diagnostic  of  transverse  pressure 
by  employing  the  same  techniques  to  measure  beam  rotation  as  used  in  HIPAC. 
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Appendix  A:  Conservation  of  Magnetic  Moment 


To  investigate  the  conservation  of  the  magnetic  moment  u  in  the  ordering 


given  in  Eqs.  (7)-(14),  we  consider  the  two  possible  effects  that  might  break 


the  invariance  of  p.  The  first  is  the  variation  of  the  relativistic  cyclotron 


frequency  u>  / y  due  to  drift  motion.  That  is,  we  must  estimate  5o>  /u>  , 

c  cr  cr 


where  5  represents  the  variation  because  of  guiding  center  drift  into  a  field 


imhomogeneity  within  one  relativistic  cyclotron  period  and  c o  ■  to  h • 

cr  c 


We  first  write  5<o  /u»  _  •  5<o  /to  -  6y/y*  Since  Ba  ■  g(Pfl)/r,  we  find 

cr  cr  c  c  o  tj 


*V  St  V 

’&  '  R  u)  R’ 

c  c 


(A-l) 


which  is  of  order  e 


(6r  is  the  change  in  radius  r  during  one  cyclotron 


period.)  Next,  using  y»c  +  q$  ■  E(PQ),  find 


(A— 2) 


where  the  first  of  the  two  factors  on  the  right  in  (A-2)  is  due  to  the 


relative  variation  of  $  along  a  drift  surface  due  to  toroidal  effects  and  the 


second  factor  is  the  angle  swept  out  along  the  drift  surface  in  one  cyclotron 


period.  From  q$/Ymc^  -  1/y  ->  we  see  that  (Sy/y)^  is  of  order  z* . 


The  other  factor  which  might  influence  invariance  of  the  magnetic  moment 


is  variation  of  the  relativistic  cyclotron  frequency  ojc/y  over  one  gyroradius 


across  the  drift  surfaces.  Again  we  have  two  contributing  factors,  the  first 


of  which  is 


6B0  WPg)]  $r 
"c  «  '  B9  "  "  R  ’ 


(A-3) 
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*  •  *  .  1 
•V-* 


1 


*.  *w 


The  first  term  is  due  to  the  diamagnetic  depression  from  the  poloidal  beam 

currents,  and  the  second  term  is  due  to  the  overall  1/r  dependence  of  B^, 

i.e.  BQ  -  g/r.  We  find  6g/g  (Ag/g)(p/a),  where  Ag/g  is  the  relative 

2*4"  2  6 

depression  in  self  across  the  whole  beam,  and  is  of  order  e 

Therefore,  we  find  5g/g  -  e®/3+25.  Also,  6r/R  (p/a)(a/R)  -*  Finally, 

we  obtain 


(ii)  .  p , 

{y  ;g  2  a’ 

'fmc 


(A-4) 


where  the  first  factor  is  the  drop  in  electrostatic  potential  from  the  center 

1+6 

to  the  edge  of  the  beam,  of  order  1/y  -  e  .  Therefore,  we  find 

.  5/3+6 

C6y/y)  -  e 
© 

Summarizing,  we  find  that  the  total  variation  6u  /u>  over  one  cyclotron 

cr  cr 

period  is  dominated  by  the  variation  in  y  across  the  drift  surface  because  of 
the  finite  Larmor  radius  of  the  electrons,  and  is  order  +  this  Is 

small  for  5  in  the  range  of  Interest  (-1/3  <  6  <  1/3)  so  that  the  magnetic 
moment  is  a  good  invariant. 
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